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Abstract 

In this paper we compute new reaction rates of the Smoluchowski equation which takes account 
of correlations. The new rate K = K MF + K c are the sum of two terms. The first term is 
the known Smoluchowski rate with the mean-field approximation. The second takes account of a 
correlation between clusters. For this purpose we introduce the average path of a cluster. We relate 
the length of this path to the reaction rate of the Smoluchowski equation. We solve the implicit 
dependence between the average path and the density of clusters. We show that this correlation 
length is the same for all clusters. Our result depends strongly on the spatial dimension d. 

The mean-field term Kfj F = (Dj + Dj)(rj + rj) d_2 , which vanishes for d = 1 and is valid up 
to logarithmic correction for d = 2, is the usual rate found with the Smoluchowski model without 
correlation (where ri is the radius and Di the diffusion constant of the cluster). We compute 
a new rate: the correlation rate Kf- = (Di + Dj)(rj + ri) d ~ 1 M(^j^) is valid for d > l(where 
M(a) = Ylt=l i a Ni is the moment of the density of clusters and df is the fractal dimension of the 
cluster). The result is valid for a large class of diffusion processes and mass-radius relations. This 
approach confirms some analytical solutions in d = 1 found with other methods. We show also 
Monte-Carlo simulations which illustrate some exact new solvable models. 
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There are different microscopic processes which lead to the model with the Smoluchowski 
equation. A comprehensive review can be found in pjand the derivation of the reaction 
rate can be found in with the mean-field approach. Particles can be injected into the 
system during time evolution [ 3 * 4 * 5 * 6 * 7 ], the injection of particles can be periodic in time 
P^. For a very anisotropic surface the model is practically one-dimensional [ 10 | 11 | |. Hence 
we will assume the most general model for an arbitrary space dimension d and flux F(t) of 
particles injected into the system. Our model has two properties: 



1. Brownian diffusion: clusters of mass k diffuse with a Brownian motion with a diffusion 
constant 

2. irreversible aggregation: clusters interact through a contact process, namely two par- 
ticles in contact aggregate irreversibly. 

The Smoluchowski approach disregards the microscopic density rik(t,r) of a cluster of 
mass k and considers the macroscopic density, i.e. space average density: Nk(t) = 
j ^ d ) r J d^rn(t, r). The Smoluchowski equation for the density of a cluster is: 

J AT -i +0° 

= j E KijWi ~ N * E ***** + F (*)*M. (!) 

i+j=k i=l 

where is the reaction rate and 5jj the Kronecker symbol. We assume the following 
initial conditions: if there is no source F(t) = we have N^t = 0) = Si^, otherwise with a 
source we have Ni(t = 0) = 0. A useful quantity is the moment of order q: 

The time evolution of the moment is: 

M g (t) = \ J2i tj K i;j ((i + jY -i*- f)N t N 3 + F(t). 

For q = 1 the moment is the total coverage Mi, i.e. the total mass, and for q = it is the total 
cluster density M . Without source the total mass is conserved M 1 (t) = Mi(0) = N^O) = 1. 
The calculation of the reaction rate depends on some physical assumptions The mean-field 
theory assumes we can exchange the many-body problem by a one-body problem with a 
convenient external field. For this purpose we define asymmetric rates, the rate that a cluster 
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of mass i aggregates with the cluster j: K^j (the reaction rate is a symmetric function of 
K^j). Following p], (we chose a characteristic length instead of a characteristic life- 
time used in p}) and dimension analysis, we define the average length of a cluster with the 
following equation: 
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the life-time 7$ of the cluster is defined with the relation £j = \fDjji p). We follow the 
mean-field approach p); we chose an arbitrary particle n»(r, t) and the microscopic equation 
of the cluster reads 

X ~ i>: A%„ (3) 



where A'^ is the symmetrical Laplace operator in d dimensions: 

A(«0 



d-i a 

r dr ' 



The fluctuations (JVj — n$) around the average value iV, are proportional to the variation 
of the microscopic density A^n^. We consider a spherically symmetric solution of Eq. 
The external field is the same as for the mean-field approach [ 2 ] , the boundary conditions 
rii(r — > +00) = iVj and n(r = r^) = express the convenient external field. Physically it 
means that close to a cluster the microscopic density vanishes, i.e we have a perfect sink 
and far away from a sink the microscopic density is close to the average density (we assume 
an early stage of clustering). The stationary spherically symmetric solution of the Eq. ^ ( 
we chose as origin the contact between the cluster i and the point-like cluster j) reads |] 
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where K is the modified Bessel function p=). The reaction rates is proportional to the flux 

Didni\ 
J AT; Or l r=r j- 



of cluster i into a cluster j [ 2 J J". drj^^-| r=ri ., then: 
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We have used the asymptotic expansion of the modified Bessel Function ~ 1 + 2 for 

x — ■> +00 [ 14 J( on length scale of the order <C Jj-, i.e. x — > we have a logarithmic 
behavior ~ — which is typical at the critical dimension d — 2). 

The rate for all dimensions is formally: 

^ * 2 • 1 )• ( 4 ) 

Si 

Eq. (4) contains two terms. The first term ?"j~ 2 , which is (strictly speaking) valid only for 
d > 3, we call the mean-field term. The second term which is valid only for d > 2, contains 
the correlation, implicitly defined in £j, we call the correlation term. Eq. Q is also valid for 
d = 1 if we put Tj~ 2 = 0. One then combines the Eq. (4) with Eq. (2j), we get an equation 
of the second order on 

j i 

The approximate solution, for a large time scale where the correlation length is short, i.e 
ii < r jt is 

^ = 

Notice that the average length £jis independent of the mass i. It is on average the inverse 
of an effective coverage 

Erf 1 ** 

and the life-time r% = ^ depends on the mass i. We obtain the asymmetric rate: 
The symmetrization of the reaction rate leads to: 

+00 

K i}j = (A + DMir, + n) d - 2 + (r, + r^^rf" 1 ^). (6) 

fc=i 

One generally assumes that the gyration radius scales with the mass ri = i a ( in appropriate 
length units), where a = 1/df is the inverse of the fractal dimension of the boundary of the 
cluster( for compact island df — d — 1, and df = 1.72 for Diffusion Limited Aggregation 
process in in d = 2 | 12 j). The effective coverage is for compact islands, M^-i)/dt * s close to 



the coverage M\ when d 3> 1. 

From the Eq. ([6| we define two rates: the mean-field rate K^ F = (A + Dj)(rj + r^) d ~ 2 and 
the rate with correlation Kf - = (A + Dj)(rj + r i ) d ~ 1 M a ( d _ 1 ) . The difference between these 
both rates is that we exchange a spatial dimension (rj + r$) in mean-field rate through a 
moment M a (d-i) in the rate with correlation. We have hence a Taylor expansion around the 
correlation number M "l d ~ 1) . In d = 1 only the correlation term remains and the rate is: 

Kij = (A + ^i)M . (7) 

We assume M a ^-\) — > M when d — > 1, because the scaling effect of the radius of gyration 
vanishes in d = 1, i.e. the boundary for an aggregating particle is point-like. 
For a time scale larger than t c , defined by M a (d-\) (t c ) ~ ( r i + r i) ; the correlations dominate, 
hence the reaction rate is essentially: 

= Kg = (A + ^Ofo + r^M^y (8) 

We can derive solutions of the Smoluchowski equation with only rate with correlation 
Eq. ^ from known solutions with mean-field rate Kf lF = (A + Dj)(rj + ri) d ~ 2 in a 
special case, without source (F = 0) and with a rate given from Eq. ([8])the solution of the 
Smoluchowski equation in d dimensions is equivalent to the solution with the mean-field 
rate K MF without source (F = 0), up to a time rescaling dt = M a u-i)(t)dt, in d + 1 
dimensions. 

With our method a large class of processes in d = 1 can be solved analytically. 
In order to confirm this approach we present some analytical solutions and numerical evi- 
dence. 

1. The solution of the Smoluchowski equation with K it j = 1 and without source (F = 0) is 
known analytically p], ~ t~ 2 and M ~ t" 1 , hence the time rescaling dt = M (t)dt 
leads us to the known solution in d = 1 of the clustering process with mass-independent 
diffusion. This solution was found with another method p*) is Mq ~ t~ l l 2 . Notice 
that through our method we obtain the asymptotic behavior for all clusters, Nk ~ t^ 1 
( F1GQ which was not possible with the method used in p^|. 

2. With constant source F = const, and = 1 we have the following equations: 
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FIG. 1: Cluster density N± as a function of time (for comparison the red curve is ~ t 1 ) 




FIG. 2: Asymptotic value of the total cluster density Mq as a function of the flux F (for comparison 
the red line has a slope F 1 / 3 ) 

Nk = M -£ i+j=k NiNj - 2N k M 2 + 5 1>k F and M = -M 3 + F. 

The asymptotic solutions are M ~ F§ ( FIG (2} and iV fc ~ fc _ § ( FIgQ . For 
comparison the solution with the mean- field term pSJin d = 2 is JV ~ Ft and 

3. For the sake of completeness we present an analytical solution already known: D k = 
and source F = const. The solutions were found with another method in 



FIG. 3: Asymptotic value of the cluster density Nk as function of the cluster mass /c(for comparison 
the straight line ~ /c~ 3 / 2 ) 

are Ni ~ t~ 1 / 2 and total cluster density Mq ~ t 1 / 4 . With our method we obtain the 
same rates and therefore the same solution. For comparison the mean-field approach 
gives N\ ~ t^ 1 / 3 and Mq ~ i 1 / 3 which are different to the numerical solution and the 
experimental result. 

Hence we have shown that this approach leads to the correct description of irreversible 
aggregation with brownian diffusion in one dimension for different models. We hope to find 
similar results in d > 1 beyond the time scale where the mean-field rate is valid, i.e. for 
t » t c . 

There are some open questions: 

1. Higher order of the Taylor expansion: can we compute other rates when the correlation 
between cluster is more important, i.e. can we obtain the full Taylor expansion around 
the correlation parameter M ° ,< - d - 1) Perhaps larger correlation effect can be obtained by 
modifying the boundary condition of Eq. ^ 

2. multiple particles reaction: how can this approach be generalized to systems where 
the binding energy is so tight that binary reaction is negligible and we have to consider 
the reaction A\ + ■ • ■ + A n — > B\ H + B m with an arbitrary n and m ? 
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